A new mechanism for granular segregation 
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A novel process is described that produces horizontal size segregation in a vertically vibrated layer 
of granular material. The behavior is a consequence of two distinct phenomena that are unique to 
excited granular media: vibration which causes the large particles to rise to the top of the layer, 
and a vibrating base with a sawtooth surface profile which can produce stratified flows in opposite 
directions at different heights within the layer. The result of combining these effects is that large 
and small particles are horizontally driven in opposite directions. The observations reported here 
are based on computer simulations of granular models in two and three dimensions. 

PACS numbers: 45.70.Mg, 02.70.Ns 



Granular media have been said to embody properties 
of solids, liquids and gases, but in actual fact studies of 
the granular state continue to produce surprising results 
that are unique to this class of matter jl], |^, |) . Computer 
modeling plays an important role in these developments, 
due to the fact that once it has demonstrated an ability to 
reproduce the observed behavior it can be used in trying 
to understand the mechanisms involved; examples of the 
dynamics of granular media that have been examined 
by this kind of approach include vertical size segregation 
by both vibration || and shear |6[ f?J, and vibration- 
induced surface waves [g, |l(| |ll| . 

A further instance of unusual behavior in granular mat- 
ter that has come to light very recently involves a granu- 
lar layer that is vibrated vertically by a base whose sur- 
face profile consists of sawtooth-like grooves (in the ex- 
perimental realization the material is contained in a nar- 
row annular region between two upright cylinders and 
is thus practically two-dimensional). What is observed 
fl2|] , both experimentally and in computer simulation, is 
the occurrence of horizontal flow. Moreover, the flow di- 
rection and magnitude depend on many of the param- 
eters required to specify the system in an apparently 
complex manner. Subsequent more detailed simulation 
of this phenomenon [|l3| revealed that the induced flow 
rate actually varies with height within the granular layer; 
even more surprising is the observation that oppositely 
directed flows can exist simultaneously at different lev- 
els in the layer. Once the strongly stratified nature of 
the horizontal flow is appreciated, the sensitive param- 
eter dependence of the overall horizontal motion of the 
layer is readily understood as being a consequence of the 
competition between opposing stratified flows. 

A combination of two of the phenomena that are 
unique to granular matter, namely the already famil- 
iar vertical size segregation under vibration, the so-called 
"Brazil nut" effect, and the recently observed stratified 
flow outlined above, suggests an entirely new mechanism 
for separating the components of a granular mixture ac- 



cording to particle size. If, under suitable vibration con- 
ditions, the upper and lower layers of the material are 
able to move horizontally in opposite directions under 
the influence of a sawtooth-shaped base, and the same 
vibrations cause the large particles to climb towards the 
top of the layer, then it is apparent that some degree of 
horizontal segregation of large and small particles should 
occur as a consequence. A process based on behavior 
of this kind could well be significant from an industrial 
perspective. In order to determine whether this pro- 
posed segregation actually occurs, computer simulations 
based on a discrete-particle model for granular material 
have been carried out in both two and three dimensions. 
Most of the results are for the former, which is essentially 
equivalent to a very narrow upright container, since con- 
siderably smaller numbers of particles are involved in the 
calculations. 

The granular model employed here is the same as in 
the horizontal flow simulations p3| ; early work with sim- 
plified, inelastic, soft-particle models of this kind ap- 
peared in j^, Q , and a number of similar models are now 
in widespread use for modeling the dynamics of gran- 
ular flow. The interaction used to provide the spheri- 
cal shape of the granular particles has a Lennard- Jones 
(LJ) form, truncated at a range r c where the repul- 
sive force falls to zero. For particles located at n and 



Tj this is flj 

for m < r c ■■ 
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(<7j + <Jj)/2. <7i is the approximate diameter of particle 
i, although because the particles are slightly soft the di- 
ameter is not precisely defined. Normal and transverse 
damping forces (the latter is optional) act during the col- 
lision (i.e., while r, 
-Jn{r i3 ■ 
fii = -min^l/, 



< r c ); the normal damping force 
''fa, and the transverse force is 
"z)?„ , where is the rela- 



*Electronic address: rapaport@mail.biu.aci 



is fij — mv »j " 1 *3 >' *j ' 

V = -min(/z|/r.+/^|,7 s |Uy l/ „^, „ y 
tive transverse velocity at the point of contact (which de- 
pends on the angular velocities of the particles and their 
relative translational velocity). The value of the static 
friction coefficient is /i = 0.5, and the normal and trans- 
verse damping coefficients are j n = 7s = 5. Particles are 
also subject to a gravitational acceleration g. 

For convenience the results are expressed in terms of re- 
duced units, in which length measurements are based on 
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the diameter of the small particles. Such units are readily 
converted to actual physical units. If, for example, the 
diameter is actually I0~ 3 m, then the reduced time unit 
consistent with the value g = 5 used subsequently corre- 
sponds to 0.022s, so that the vibration frequency value 
/ = 0.4 employed in these simulations is equivalent to 
rs 18Hz, typical of the values used experimentally. Hor- 
izontal flow rates |fl3f within the strata are then of order 
10 _2 m/s. 

The particle size distribution is bimodal, with values 
randomly distributed over narrow ranges extending from 
the nominal small and large sizes to values 10% lower. All 
particles have the same material density. The sawtooth 
base for the two-dimensional simulations is constructed 
(for reasons of computational convenience) from a set of 
grainlike particles positioned in a manner that describes 
the desired asymmetric sawtooth profile; the base parti- 
cles are forced to oscillate vertically in unison to produce 
the effect of a sinusoidally vibrated base, and it is these 
vibrations that drive the system. The base particles in- 
teract with the freely moving granular particles using the 
same force laws as above; in order to ensure reasonably 
straight sawtooth edges the diameter of these particles is 
0.33 and the distance between their centers is 0.17. Fur- 
ther discussion of the model appears in and other 
more general aspects of the computational methodology 
are to be found in Jl5| . Unlike the earlier work, which 
dealt with a horizontally periodic system, particles here 
the are confined to a container with reflecting side walls; 
the container is sufficiently high that particles never reach 
the upper boundary. 

Figure [l] shows a series of screen snapshots (the image 
resolution is limited) for a mixture in which there is a 
15% concentration of large particles with diameter 1.4. 
The first image corresponds to an early state of the sys- 
tem, with large particles randomly located throughout 
the layer. Subsequent images show the particle distri- 
butions after 1000, 2000, 4000 and 16000 vibration cy- 
cles; the occurrence of size segregation, in this instance 
marked by the leftward motion of the large particles, is 
clearly visible. 

The presence of the side walls interferes with the strat- 
ified horizontal flows observed with periodic boundaries 
(corresponding to an effectively infinite system) |l3| ] ; nev- 
ertheless, in a system of the width considered here, the 
stratified flow not too close to the walls is sufficiently 
strong to drive the segregation. Any practical implemen- 
tation of such a separation mechanism would of course 
require some means of continuously injecting a granular 
mixture at a point above the layer somewhere near its 
center, and extracting the (at least partially) segregated 
products from opposite ends of the container through 
suitably located openings. 

Concentration profiles for the large and small parti- 
cles at various times are shown in Figure The gradual 
emergence of segregation is apparent, beginning with an 
initially uniform system, and finishing with practically 
all the large particles having migrated to the left side 




FIG. 1: Screen images showing stages in the size segregation 
process for a two-dimensional system; the large particles are 
darkly shaded. 
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FIG. 2: Variation of the local concentration of large and small 
particles (solid and dashed curves respectively) across the sys- 
tem at several stages during the segregation process. 



of the container. The system width is 180 (in reduced 
units) for the particular case shown here, and the nomi- 
nal number of layers is 8. The frequency is / = 0.4 and 
the amplitude A = 1. The dimensionless acceleration 
r = (2-7T f) 2 A/ ' g is an important quantity in vibrating 
systems, with r ~ 1 the minimum required to excite the 
layer; a value T = 1.26 is used here, obtained by setting 
g = 5. The base contains 20 sawteeth of height 2, with 
an asymmetry such that the right edge of each tooth is 
practically vertical. For the case of granular particles 
with a unimodal size distribution fill this choice of pa- 
rameters produced oppositely directed stratified flows of 
similar magnitude (there the system width was just 90, 
so the relevant velocity profile is for 10 sawteeth). The 
large particle size is 1.2 and the concentration 15%. Re- 
sults are averaged over short time intervals and over six 
independent runs (with different initial states), and are 
smoothed to reduce fluctuations. Similar behavior is ob- 
served if the LJ overlap potential used here is replaced 
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FIG. 4: Trajectories of 2% of the particles after segregation 
has ended (large particle trajectories are in black), for pairs 
of values of large particle size and number of sawteeth (1.4, 
20), (1.2, 20), and (1.2, 80). 



FIG. 3: Position of the center of mass of the large particles as 
a function of time for different particle sizes, and for 20 and 
40 sawteeth (solid and dashed curves). 



by linear or Hertzian potentials; thus the results do not 
depend on the choice of interaction. The fact that segre- 
gation occurs for a large particle size of only 1.2 suggests 
a highly sensitive mechanism. 

In the case of the unimodal size distribution, the strat- 
ified flows were observed to depend, often in a com- 
plex manner, on the various parameters used to specify 
the system. A similar dependence also occurs for mix- 
tures. Space does not permit a detailed analysis of the 
dynamical phase diagram, which would show how the 
separation rate and direction (as well as the purity of the 
segregation products) depend on each of these param- 
eters. Instead, a few selected examples illustrating the 
variation with certain key parameters will be presented, 
with emphasis on cases where segregation is strong. 

Particularly important from a practical perspective is 
the efficiency of the segregation mechanism for different 
large particle sizes. Figure |^ shows this size-dependence 
for two sawtooth widths. The degree of segregation is 
expressed in terms of the time-dependent position of the 
center of mass of the large particles; the occasional par- 
ticle trapped in a sawtooth groove before its migration 
is complete will adversely affect such a measurement, as 
will any other effects preventing the formation of a well- 
defined vertical boundary between the two segregated 
species (see below). The results are again averaged over 
six independent runs. The final value depends on the 
large particle size; in the case of 20 sawteeth it increases 
monotonically with size, reflecting the increasing space 
occupied by the large particles, but for 40 sawteeth the 
results for diameter 1.2 deviate from the expected se- 
quence. The times required to reach the limiting position 
are essentially the same. 

While a detailed picture of the particle movement over 
the entire history of the segregation is not readily ob- 
tained because of the strong degree of mixing that occurs, 
a study of the trajectories for a subset of particles, once 
the process has effectively ended, demonstrates where the 



particles tend to be located. Figure |4| shows short tra- 
jectory segments for just 2% of the particles after 35000 
vibration cycles; each trajectory consists of five points 
at 100 cycle intervals. In the first two examples (with 
20 relatively wide sawteeth, and large particle sizes 1.4 
and 1.2) the large particles congregate on the left; they 
occupy a region that includes the section of the layer 
abutting the wall extending from top to bottom, and ad- 
jacent to this a roughly wedge-shaped region that tapers 
upwards to the right. The non- rectangular shape of this 
region explains why the center of mass position does not 
fully characterize the final state. In the third example 
(80 very narrow sawteeth, particle size 1.2) segregation 
occurs in the opposite direction and is less complete. 

In those cases where large particles migrate to the left, 
the shape is a consequence of a dynamic equilibrium in- 
volving both stratified flow, in which particles in the 
upper levels are driven horizontally to the left by the 
sawteeth and lower particles to the right, and Brazil-nut 
"buoyancy" which raises the large particles. In a practi- 
cal implementation of such a segregation apparatus, large 
particles would be removed near the upper left corner and 
small particle via an opening near the lower right. The 
case where large particles migrate to the right reflects the 
absence of leftward stratified flow due to the sawtooth 
base. In all cases shown here, irrespective of direction, 
the horizontally driven large particles eventually displace 
the small particles from the appropriate side wall. 

The corresponding three-dimensional version of this 
system consists of a rectangular base with linear 
sawtooth-like grooves; similar results would be expected 
from simulations of such a system. Since it is more in- 
teresting to consider a problem involving an alternative 
geometry, a system with circular symmetry is described 
here. 

In this three-dimensional example the base consists of 
a set of concentric circular grooves, shaped to produce 
a radial profile identical to that used in two dimensions. 
Interactions between particles and the base are handled 
with a simple extension of the approach used above, by 
employing a series of concentric donut-like base "parti- 
cles" to represent the grooves (only the radial coordi- 
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FIG. 5: Segregation in a three-dimensional system; a wedge- 
shaped region has been removed from the front and the con- 
tainer is not shown (large particles are darker). 

nates of such particles enter into the force computation). 
Since a substantially larger number of granular particles 
are necessary to achieve a linear size similar to that of 
the two-dimensional system, only a limited exploration 
of this problem has been undertaken. The transverse 
damping force is omitted here; while its absence does not 
seriously influence the observed behavior, indicating once 
again that the phenomenon is not especially sensitive to 
the details of the model, it does reduce the computational 
effort. 

Figure is a screen image from one such simulation 
showing the system after 2000 vibration cycles. The con- 
tainer diameter used here is 180 and the nominal layer 



thickness is 8; there are 11 concentric sawtooth grooves, 
and the large particles have size 1.5 and concentration 
33%. The sawtooth profile, whose influence is apparent in 
the positioning of the particles closest to the base, is ori- 
ented so that large particles should migrate towards the 
outer boundary of the cylindrical container while small 
particles move inwards, which indeed is exactly what is 
seen to occur. 

In summary, computer simulations employing a type 
of granular model whose viability has been established 
in other studies of grain flow have been used to inves- 
tigate a novel mechanism for achieving size segregation. 
Selected results have been presented here, with a more 
detailed determination of the parameter dependence to 
be published in due course. Since separation and sorting 
are essential functions in the processing of bulk granular 
materials, such an approach, implemented for example as 
a continuous flow device, has potential industrial value. 
It remains to be seen whether the behavior of real gran- 
ular matter is in accord with the predictions of these 
simulations. 
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